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Abstract 

In this review, we discuss recent experiments that investigate how the quantum sate of a superconducting qubit 
evolves during measurement. We provide a pedagogical overview of the measurement process, when the qubit is 
dispersively coupled to a microwave frequency cavity, and the qubit state is encoded in the phase of a microwave tone 
that probes the cavity. A continuous measurement record is used to reconstruct the individual quantum trajectories 
of the qubit state, and quantum state tomography is performed to verify that the state has been tracked accurately. 
Furthermore, we discuss ensembles of trajectories, time-symmetric evolution, two-qubit trajectories, and potential 
applications in measurement-based quantum error correction. 
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1. Introduction 

The standard description of quantum mechanics con¬ 
siders the time-evolution of isolated quantum sys¬ 
tems whose unitary dynamics are governed by the 
Schrodinger equation. Measurement is treated as an in¬ 
stantaneous non-unitary process through which a quan¬ 
tum system is projected into an eigenstate of the mea¬ 
sured observable with a probability given by Bom’s 
rule. In reality, no system is completely isolated from 
its environment, and measurements are never traly in¬ 
stantaneous, but occur over some finite timescale de¬ 
termined by the details of the interaction between the 
measured system and its environment. The theory of 
quantum trajectories inia considers measurement as a 
continuous process in time, describing how the state of 
the quantum system evolves during measurement. 

Due to the intrinsic quantum fluctuations of the envi¬ 
ronment, measurement is an inherently stochastic pro¬ 
cess. If a quantum system starts in a known quantum 
state |i/'(0)), then by accurately monitoring the fluctua¬ 
tions of its environment it is possible to reconstruct sin¬ 
gle quantum trajectories which describe the evo¬ 
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lution of the quantum state in an individual experimen¬ 
tal iteration. 

The concept of quantum trajectories was first devel¬ 
oped in the early 1990’s as a theoretical tool to model 
continuously monitored quantum emitters iHlllll. For 
the next decade, quantum trajectories were used pri¬ 
marily in the quantum optics community, as a theoret¬ 
ical tool for numerical simulations of the ensemble be¬ 
havior of open quantum systems Ena. Typically, the 
master equation of an open quantum system cannot be 
solved analytically, and thus numerical solutions are of¬ 
ten necessary. For a Hilbert space of dimension N, the 
density matrix p consists of real numbers, and the 
computational time required to solve for its time evolu¬ 
tion through the master equation scales as i). In 
contrast, the pure quantum state \>f/{t)) of an individ¬ 
ual quantum trajectory can be described by N complex 
numbers. Therefore, it is often advantageous to simulate 
an ensemble of stochastic quantum trajectories, which 
can be averaged together to recover the evolution of the 
density matrix, p{t). Although the formalism of quan¬ 
tum trajectories is constructed from standard quantum 
mechanics Q, it can provide insight into foundational 
questions such as the quantum measurement problem 
El 13 [To) and bears a close resemblance to the consis¬ 
tent histories interpretation of quantum mechanics IfTTl . 
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Despite widespread theoretical use, quantum trajec¬ 
tories have only been investigated in a handful of ex¬ 
periments, due in part to the difficulty of perform¬ 
ing highly efficient continuous quantum measurements. 
The earliest experiments to continuously monitor indi¬ 
vidual quantum systems were in the regime of strong 
measurement, where the system is quickly projected 
into an eigenstate of measurement, destroying any in¬ 
formation about the phase of a coherent superposition. 
In such experiments, it is possible to track the ‘quan¬ 
tum jumps’ between eigenstates |[I21|I3][I1[T51. Cav¬ 
ity quantum electrodynamics (CQED) experiments with 
Rydberg atoms have explored the weak measurement 
regime, tracking the quantum trajectories of a cavity 
field as it collapses from a coherent state into a pho¬ 
ton number eigenstate ifThll . Other CQED experiments 
have used a cavity probe to continuously track the posi¬ 
tion of individual Cesium atoms C3. Quantum trajec¬ 
tories were first considered for solid state systems in the 
context of a quantum dot qubit monitored in real time 
by a quantum point contact charge sensor IHIIIII]. In 
2007, the conditional measurement dynamics of a quan¬ 
tum dot were investigated experimentally ll20l . More 
recently, quantum trajectory theory has been used to 
solve for the conditional evolution of a continuously 
monitored superconducting qubit IITTI l22l . These re¬ 
sults, when combined with recent advances in nearly- 
quantum-limited parametric amplihers, which can be 
used to achieve highly efficient qubit readout, have en¬ 
abled a detailed investigation of measurement backac- 
tion Enm. 

In this article, we review recent experiments ll25l \26l 
IZTl | 28 l which, by weakly probing the held of a mi¬ 
crowave frequency cavity containing a superconducting 
qubit, track the individual quantum trajectories of the 
system. These are the hrst experiments, on any system, 
which use quantum state tomography at discrete times 
along the trajectory to verify that the qubit state has been 
faithfully tracked. Erom the perspective of quantum in¬ 
formation technology, these experiments demonstrate 
the great extent to which the process of measurement 
is understood in this system, and may inform future ef¬ 
forts in measurement-based feedback Il29l [30] |3T] |32l 
for state stabilization ll^ and error correction. 

The review is organized as follows. In section 2, 
we present a physical picture for continuous quantum 
measurement of superconducting qubits. In section 3, 
we demonstrate how to use a measurement result to 
reconstruct the conditional qubit state after measure¬ 
ment. In section 4, we explain how to reconstruct 
and tomographically verify individual quantum trajec¬ 
tories. Then, in section 5, we examine ensembles of 


quantum trajectories to gain insight into qubit state dy¬ 
namics under measurement. In section 6, we discuss 
time-symmetric evolution under quantum measurement, 
and in section 7 we demonstrate quantum trajectories of 
a two-qubit system. Einally, in section 8 we explore 
potential applications in measurement-based feedback 
control and continuous quantum error correction. 

2. Continuous measurement of superconducting 

qubits 

The experiments discussed in this review use artifi¬ 
cial atoms formed from superconducting circuits. We 
focus in particular on the transmon circuit 1341 (Eig. 
lA) which is composed of the non-linear inductance of 
a Josephson junction and a parallel shunting capacitance 
Cj;. This circuit is characterized by the Josephson en¬ 
ergy scale Ej = MqIIc and the capacitive energy scale 
Ec = e^/2Cj;, where Iq is the junction critical current 
and e is the elementary charge. A typical transmon cir¬ 
cuit, with Ej/h ~ 20 GHz and Edh ~ 200 MHz, has 
several bound eigenstates (Eig. IB) with energies E,„, 
where m is a whole number that indexes the states. The 
lowest two levels form a qubit subspace, with transition 
frequency woi/27r s {Ei - Eo)lh ~ 5 GHz, and the dif¬ 
ference in frequency between transitions to successively 
higher levels is given by the anharmonicity a ^ Eq- 
Due to the large EjjEc ratio the transmon qubit is in¬ 
sensitive to charge noise, which, when combined with 
low-loss materials lT5l and designs that minimize 
the participation of surface dielectric loss El, has al¬ 
lowed for planar qubits with coherence times of many 
tens of microseconds. 

In order to to control the qubit’s interaction with its 
external environment, it is coupled to the fundamen¬ 
tal mode of a three dimensional waveguide cavity of 
frequency u)c at a rate g, realizing a cavity quantum 
electrodynamics (CQED) architecture. In the dispersive 
regime, where the qubit-cavity detuning coq - a>r is large 
compared to g, the system is described by the Hamilto¬ 
nian H - EIq + //int -I- //r, where 

Hmi = -fixa^acr^ (1) 

77r = ^^o-y. (2) 

Here Hq describes the uncoupled qubit and cavity en¬ 
ergies and decay terms, h. is the reduced Plank’s con¬ 
stant, X is the dispersive coupling rate, d and a are 
the creation and annihilation operators for the cavity 
mode, and cr^ is the qubit Pauli operator that acts on 
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Figure 1: Dispersive measurements of a superconducting qubit. (A) A transmon qubit couples dispersively to a microwave frequency cavity. 
Signals that reflect off of the cavity are amplified by a nearly quantum-limited lumped-element Josephson parametric amplifier. (B) Schematic 
representation of the transmon potential and corresponding energy levels. The two lowest energy levels form the qubit subspace. (C) Reflected 
signal as a function of probe frequency. The resonance frequency is shifted by 2x depending on whether the qubit is prepared in the |0) or 11) state. 
{D,E) The cavity is probed with a coherent microwave tone at a frequency cOm, initially aligned along the Xi quadrature. After leaving the cavity, 
the tone acquires a quit-state-dependent phase shift. (F) Phase-sensitive amplification along the X 2 (top) and the Xi (bottom) quadratures. 


the qubit state in its energy eigenbasis. Hint is an inter¬ 
action term which equivalently describes a qubit-state- 
dependent frequency shift of the cavity of -xo'z (with 
the |0) state dehned as cr^ = h- 1) and a qubit frequency 
that depends on the intracavity photon number « = a^'a 
(an a.c. Stark shift). describes the effect of an op¬ 
tional microwave drive at the qubit frequency which 
causes the qubit state to rotate about the y axis of the 
Bloch sphere at the Rabi frequency Q. 

Because the Hint commutes with cr^, the qubit-state- 
dependent phase shift can be used to perform continu¬ 
ous quantum non-demolition (QND) measurment of the 
qubit state in its energy eigenbasis 13^ [39)1 . Figure 1C 
illustrates the phase of the reflected signal as a function 
of frequency. If we choose to measure at a frequency 
u),„ = (w|o> + W|i))/2, where W|o> and W|i) are the cavity 
frequencies when the qubit is in the ground and excited 
states, respectively, then the phase difference in the in¬ 
ternal cavity field for the two qubit states is given by 
A6 = 4[;y|//c, where k is the cavity decay rate. In the 
experiments presented here, we work in the small phase 
shift limit, with \x\Ik ~ 0.05. 

We probe the cavity by applying a measurement tone 
at frequency a),„ initially aligned along the Xi quadra¬ 
ture (Fig. ID). Due to the vacuum fluctuations of the 
electromagnetic held, the quadrature amplitudes Xi and 
X 2 of this held will Huctuate in time. The circle in Fig¬ 


ure ID represents the Gaussian variance of the input sig¬ 
nal time-averaged for a time At. The area of the circle is 
inversely proportional to At. After rehecting off of the 
cavity, the measurement signal acquires a qubit-state- 
dependent phase shift, as depicted in Figure IE. In the 
small \x\Ik limit, the X 2 quadrature of the rehected sig¬ 
nal signal contains information about the cavity phase, 
which is proportional to the qubit state. Likewise, the 
X\ quadrature contains information about the amplitude 
of the cavity held and thus the Huctuating intracavity 
photon number. 

In order to track the qubit state through an individ¬ 
ual measurement, we need to accurately monitor the 
quantum huctuations of the measurement signal, which 
are typically much smaller than the thermal Huctuations 
of the room temperature electronics that are needed to 
record the measurement result. Therefore, we must hrst 
amplify the signal above this noise Hoor. State-of-the- 
art commercial low-noise amplihers, which are based 
on high electron mobility transistors (HEMTs) and can 
be operated at 4 K, add tens of photons of noise to the 
measurement signal. Therefore, a more sensitive pre- 
ampliher is needed in order to overcome the added noise 
of the HEMT ampliher. 

Over the past few years, Josephson junction based su¬ 
perconducting parametric amplihers have emerged as 
an effective tool for realizing nearly-quantum-limited 
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amplification. Phase-preserving amplifiers such as 
the Josephson parametric converter HOl amplify both 
quadrature amplitudes evenly by a factor of Vg, where 
G is the power gain of the phase-preserving amplifier, 
and add at least a half photon of noise BTI to the signal. 
Here, we focus instead on phase-sensitive amplification 
from a lumped-element Josephson parametric amplifier 
El, where one quadrature is amplified by a factor of 
2 Vg and the other quadrature is de-amplified by the 
same factor. 

When we apply a coherent measurement tone, char¬ 
acterized by an average intracavity photon number h, its 
quantum fluctuations will cause the phase coherence of 
a qubit superposition state to decay at the ensemble de¬ 
phasing rate F = S^^h/a-. The ensemble dephasing rate 
will be the same regardless of how we choose to pro¬ 
cess the measurement signal after it leaves the cavity. 
However, the backaction of an individual measurement 
will depend significantly on our choice of amplification 
scheme. 

As depicted in Figure IF, after a measurement tone 
initially aligned along the quadrature reflects off of 
the cavity and acquires a qubit-state-dependent phase 
shift and is then displaced to the origin of the X 1 -X 2 
plane by a coherent tone, we can choose to either 
amplify the X 2 quadrature (top panel) which contains 
qubit-state information or the Xi quadrature (bottom 
panel) which contains information about the fluctuating 
intracavity photon number. Consider a qubit initially 
prepared in an equal superposition of cr^ eigenstates, say 
cr^ = +1. If we perform ideal phase-preserving am¬ 
plification of Xi, we also de-amplify the photon num¬ 
ber information, and the measurement backaction drives 
the qubit state along a meridian of the Bloch sphere 
toward one of the poles. We refer to this case as a 
z-measurement, because we acquire information about 
the qubit state in the cr^ basis. If instead we amplify 
X 2 , we also de-amplify the qubit-state information, and 
measurement backaction drives the qubit state along the 
equator of the Bloch sphere. We refer to this case as 
a 0-measurement, because we can track the phase of a 
qubit superposition state over the course of an individual 
measurement. For phase-preserving amplification both 
types of backaction are present. 

We first focus on the case of the z-measurement. Af¬ 
ter the measurement tone leaves the parametric ampli¬ 
fier and passes through further stages of amplification 
we demodulate the signal and record the X 2 quadra¬ 
ture amplitude as a digitizer voltage V(t). For a mea¬ 
surement of duration At, the measurement outcome Vm 
is given by the time-average of V(t). Depending on 
whether the qubit is initially prepared in the ground of 


the excited state, the Gaussian distribution describing 
the probability attaining a particular measurement out¬ 
come will be shifted by a voltage AV oc AO. In this 
review, we define a dimensionless measurement out¬ 
come r - 2VmlAV, such that the ground and excited 
state distributions are centered about r = +1, respec¬ 
tively, as illustrated in Figure 2A,D. From these dis¬ 
tributions, we define the dimensionless measurement 
strength S = {Ijaf', where a is the standard devi¬ 
ation of the dimensionless measurement distributions, 
which scales as (Af)”'''^. We also define the character¬ 
istic timescale t over which the qubit state is projected 
as the amount of time required for the measurement his¬ 
tograms to be separated by twice their standard devia¬ 
tion, T = AAtjS. When S is large (Af » t), the ground 
and excited state histograms are well separated (Fig. 
2A), and it is possible to determine the qubit state with 
high fidelity in an individual measurement. The mea¬ 
surement projects the qubit into an energy eigenstate, 
where it will remain after measurement. Instead, if S is 
small (Af < t), the ground and excited state histograms 
overlap (Fig. 2D), and an individual measurement only 
partially projects the qubit state. 

Formally, general quantum measurements (partial 
or projective) are described by the theory of posi¬ 
tive operator-valued measures (POVM) which yield the 
probability P{m) - Tr(Q,„ p Qj„) for obtaining an out¬ 
come m, and the associated back action on the quan¬ 
tum state, p —> D.„ipPl]„IP{m}, where the operators 
obey 01,0;,, = /. For example, for the projec¬ 
tive measurements - (I ± cr^)l2, the probability 
of a measurement yielding the qubit in the -i-z state is 
P(+z) — Tr(Q+ 5 ^Q+j;) = (1 H- {crj))/2. The partial mea¬ 
surements discussed in this review are described by the 

POVM (61 m, 

a = (3) 

where, l/4fl^ = Af/Tr. The cr^ term in Q.r causes 
the back action on the qubit degree of freedom, p —> 
PlrpPl-l, due to the readout of the measurement result r, 
resulting in the measurement dynamics discussed in this 
review. 

Our ability to reconstruct the qubit state after an in¬ 
dividual partial measurement is determined by the mea¬ 
surement quantum efficiency rj^. We have established 
that the noisy measurement tone contains information 
about the qubit state and cause an ensemble dephasing 
rate of F. In general, only a fraction, of this infor¬ 
mation is experimentally accessible, and the remainder 
is lost to environmental degrees of freedom. The mea- 
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Figure 2: Continuous quantum measurement. Illustrative measurement histograms for a single time-step (A,D) along with simulated measurement 
records (B,E) and con'esponding quantum trajectories (C,F). In the top panels. At = 200 ns, r = 50 ns, and n/2;r = 8 MHz, illustrating the quantum 
jumps regime. In the bottom panels. At = 20 ns, r = 150 ns, and H = 0, illustrating the diffusive regime. 


surement efficiency can be reduced from its ideal value 
of 77 ,„ = 1 by losses between the cavity and the para¬ 
metric amplifier, described by the collection efficiency 
Tjcoi and by added noise in the amplification chain, de¬ 
scribed by an the amplification efficiency rjamp- In the 
experiments presented here, !]„ - rjcoi'namp ~ 0.4. The 
measurement strength depends linearly on 77^, and for 
dispersive measurements in the small phase shift limit 
is given by 5 = 6Ax^hrim^t j k. 

We now turn our attention to continuous quantum 
measurement. In Figure 2 we illustrate quantum tra¬ 
jectories in the limiting cases of strong and weak mea¬ 
surement. We consider a sequence of n measurements 
occurring at times {f^ = kAf) for k - 0,1,. ..,77 - 1, 
which result in a set of dimensionless measurement re¬ 
sults [rk] - {ro, ri,..., r„_i). If the qubit is prepared 
in a known initial state, then we can use the measure¬ 
ment results to track the qubit state as it evolves under 
measurement, computing the set of conditional qubit 
states [qk] - {^ 0 . ■■■. ?n-i) corresponding to an indi¬ 

vidual measurement record [rk]- Here the Bloch vector 
q - {x, y, z) describes a general mixed single-qubit state 
in terms of the the components x = \i\p&x\, y = tr[pd-y], 
and z = trfpcTj], where p is the qubit density matrix. 
In the limit where r < Af, each time-step constitutes 
a (nearly) projective measurement. In the absence of 
any non-measurement dynamics, after the first time- 
step subsequent measurements will continue to project 
the qubit into the eigenstate corresponding to the initial 
measurement result. However, any additional dynamics. 


such as energy relaxation or Rabi driving, which occur 
on a timescale faster than Af will result in discontinu¬ 
ous jumps in the measurement record corresponding to 
quantum jumps ifTSl of the qubit state (Fig. 2B,C). 

In the opposite limit, where Af « t, each measure¬ 
ment will only slightly perturb the qubit state. However, 
by performing a sequence of repeated partial measure¬ 
ments such that f„_i » T, we can realize a projective 
measurement. In this case, the noisy detector signal 
can be used to reconstruct the diffusive trajectory of the 
qubit state as it is gradually projected toward a measure¬ 
ment eigenstate (Fig. 2E,F). 

3. Reconstructing the conditional quantum state 

In this section, we discuss in detail how to recon¬ 
struct the qubit state conditioned on an individual mea¬ 
surement outcome. One approach, taken in references 
II2211281 is to solve a stochastic master equation for 
conditional qubit state. Here, we instead describe phe¬ 
nomenological approach based on a Bayesian statistics 
Ea, which provides a particularly simply approach to 
single-qubit trajectories and was used in references 1^ 
and ESi . As recently demonstrated in reference ESI, 
for the a single qubit under weak measurement and 
weak Rabi driving, both approaches yield similar re¬ 
sults. 

Consider a qubit prepared in the initial state p{t = 0), 
which is weakly measured for a time Af, yielding a mea¬ 
surement result r. Here, we show how to apply Bayes 
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rale of conditional probabilities to update our knowl¬ 
edge of the qubit state after the measurement. We hrst 
focus on the case of a z-measurement, with Q = 0. 
From Bayes rule, we have 


P{i\r) 


P{r\mi) 

P{r) 


(4) 


where i describes the basis states {|0),|1)). Here, the 
initial probabilities for finding the qubit in the ground or 
excited states are given by P{Q) - pQQ{t-Q) and P{1) - 
Pii(f = 0). By expressing the measurement distributions 
P(r\i) explicitly and taking the ratio of the conditional 
probabilities in equation Q for both basis states, we 
find that 


pii(Af) ^ pii(0)exp[-(r-H if 
poo(Af) poo(O) exp[-(r - 1)2/2a2]' 

For a qubit initially prepared in the state q: - (1,0,0) 
we hnd that 



(6) 


where the superscript ‘z’ denotes a z-measurement. 

Note that thus far we have used a classical rule of con¬ 
ditional probabilities to determine how the qubit popu¬ 
lations evolve under measurement. Following reference 
ll22ll . we account for the qubit coherence through the 
phenomenological assumption that 


where the superscript ‘0’ denotes a (^-measurement. 
Figure 3 illustrates the conditional quantum state as 
a function r for a z-measurement (panel A) and a 
(^-measurement (panel B), with r = 600 ns and Af = 
400 ns. Note that the dephasing rate y due the unacces- 
sible part of the measurement signal is the same regard¬ 
less of our choice of amplification axis. In both cases, 
a measurement outcome of r = 0 will leave y and z un¬ 
changed, but X is reduced by a factor of Exp[-yAf]. 

A useful feature of dispersive CQED measurements 
is the ability to rapidly tune the measurement strength 
by changing the amplitude of the measurement tone. 
Therefore, it is straightforward to implement experi¬ 
mental sequences which combine partial and projective 
measurement. The sequence shown in Eigure 3c is used 
to implement conditional quantum state tomography to 
verify that we can accurately account for the backaction 
of an individual measurement. The qubit is prepared 
in the initial state (1,0,0), and weakly measured for a 
time Af. Then, we perform an optional qubit rotation 
(of njl about the y axis to reconstruct x, 7t/2 about -x 
to reconstruct y, and no pulse to reconstruct z) followed 
by a projective measurement. Eor a given measurement 
outcome r, we perform a tomographic state reconstruc¬ 
tion on the sub-ensemble of experimental iterations with 
similar measurement outcomes, in the range r±e, where 
e <Si 1. Eor superconducting qubits, this technique 
was first introduced in reference ll23l . which consid¬ 
ers the case of phase-preserving amplification. Shortly 
thereafter, this technique was demonstrated for phase- 
sensitive amplification l25l . 


x^ = (7) 

Here the first term enforces normalization and the sec¬ 
ond term reflects our imperfect knowledge of the envi¬ 
ronment and leads to qubit dephasing characterized by 
the rate y = F - 1 /2 t, where F = F -i-1 /T^ is the ensem¬ 
ble dephasing rate and ~ 20 ps is the characteristic 
timescale for extra environmental dephasing. 

For the case of a ^-measurement, z remains zero, and 
X and y are periodic in the accumulated qubit phase shift, 
and are given by E2\ 


x'^ = cos|^je-^^ 
/ = -sin|^je->'^'. 


( 8 ) 

(9) 

( 10 ) 


4. Tracking individual quantum trajectories 

Consider a qubit initially prepared in a known state 
qj, which undergoes a sequence of n partial measure¬ 
ments with outcomes {rk}, as described in section]^ In 
the limit where the duration At of each measurement 
approaches zero, the set of conditional states {qk} de¬ 
scribes a the quantum trajectory q{t). For simplicity, 
from here on we restrict our discussion to the case of a 
z-measurement. When w = 0, we can calculate the con¬ 
ditional quantum state at each time-step t^ - kAt from 
equations and 0 using only the initial state and 
the time-averaged measurement signal r — llkYilZ^rk. 
However, when Q > 0 the measurement dynamics do 
not commute with the Rabi drive, and therefore the or¬ 
der of the measurement outcomes matters, and r no 
longer contains sufficient information to reconstruct the 
quantum trajectory. Instead, if Af we can perform 
a sequential two-step state update procedure introduced 
in reference ll^ . For each time-step t^, we calculate 
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Figure 3: Reconstructing the conditional quantum state. (A,B) The 
conditional quantum state after a measurement result r for a qubit ini¬ 
tially prepared in the state (1,0,0), with At = 400 ns, r = 600 ns, and 
7 = 1.3 * lO^s^T Panel A depicts a z-measurement, and panel B de¬ 
picts a 0 measurement. (C) Experimental sequence for reconstructing 
the X component of the conditional quantum state. (D) To perform 
quantum state tomography conditioned on the measurement result 
r = 1.7, we average together the projective measurement outcomes 
for the sub-ensemble of measurement outcomes where r = 1.7 ± 6. 


qk by first applying a Bayesian update to the state 
to account for the measurement result Vk, and then by 
applying a unitary rotation to account for the the Rabi 
drive during the time At. Example quantum trajectories 
are shown in Figure 4A,B for Q/2;7r = 0 and 0.4 MHz, 
respectively, and t = 1.28yus. The corresponding en¬ 
semble average evolution is shown in panels C and D. 

While pervious experiments in other systems have 
reconstructed individual diffusive quantum trajectories 
ca, reference Ea was the first to use conditional 
quantum state tomography to verify that the trajecto¬ 
ries were reconstructed accurately. Here, we present 
a brief outline of the tomographic validation proce¬ 
dure. We perform a large number of experimental it¬ 
erations ending at different times tf, which are fol¬ 
lowed by a qubit rotation and a projective measurement. 
We use a single full-length experimental iteration (with 
tf - (n — l)Af to generate a target trajectory, denoted 
q(t) = (x(t),y(t),z(t)). Then, for each experimental se¬ 
quence of total measurement duration 4 (and a given 
orientation of tomography pulse), we compute the quan¬ 
tum trajectory q{t). We perform conditional quantum 
state tomography separately at each time 4 using the 
subset of experimental iterations with x(4) = x(4) + e 


Figure 4: Reconstructing individual quantum trajectories. Here, r = 
1.28^is, 7 = 2.7 X and n/2ff = 0 (A, C) and 0.4 MHz (B,D). 

Panels A and B depict the ensemble average evolution. Panels C and 
D display simulated individual quantum trajectories ending at = 
2/is, with the x,y, and z components depicted in blue, red, and black, 
respectively. The orange regions represent a matching window of 6 = 
0.05 at t = 0.66 jjs (C) and 1.48//S (D). Sample trajectories that end 
within the matching window are shown in other colors. 

and z(4) = 2(4) + e where e « 1, and we have assumed 
that y = 0. The orange shaded regions in panels C and D 
of Figure 4 represent matching windows at 4 = 0.66 ps 
and 1.28 jUS, respectively. Trajectories which fall within 
the matching window at 4 are used in the tomographic 
reconstruction of q(tk). 

5. Distributions of trajectories 

By tomographically reconstructing individual quan¬ 
tum trajectories, as discussed above and initially 
demonstrated in reference ||25]| . we have proven that we 
can accurately track the qubit state over the course of 
any individual measurement. In this section, we con¬ 
sider how quantum trajectory experiments are useful for 
building an intuition for how the qubit state is most 
likely to evolve under measurement. As discussed in 
reference 12 ^ . distributions of quantum trajectories of¬ 
fer a convenient qualitative tool for visualizing the inter¬ 
play between measurement dynamics and unitary evolu¬ 
tion. 

The greyscale histograms in panels A and B of Figure 
5 display the simulated distribution of quantum trajecto¬ 
ries for T = l.lSfis, y - 2.7 x 10“^and f2/27r = 0.4 
MHz. Note that due to the Rabi drive, the measurement 
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initially projects the qubit preferentially toward the ex¬ 
cited state (z = -1). At intermediate times a wide range 
of qubit states are possible, and after half a Rabi period 
the qubit is preferentially projected toward the ground 
state. In experiments with superconducting qubits, r 
and Q can be readily tuned, and distributions of quan¬ 
tum trajectories are experimentally accessible for a wide 
range of parameters. 

It is also possible to consider the conditional quantum 
dynamics of the sub-ensemble of trajectories which end 
in a particular quantum state, or post-selection. Pan¬ 
els C and D of Figure 5 display the distribution of tra¬ 
jectories which end in the final state xp = 0.1 ± .08, 
Zf - 0.55 ± .08. By analyzing the statistical properties 
of such distributions, it is possible to answer questions 
of broad interest in the field of quantum control. The 
experiments of reference ES\ focus on one such ques¬ 
tion; what is the most probable path through quantum 
state space connecting an initial state |i/r,) and a final 
state Itf/f) in a given time T1 

One straightforward theoretical approach to this 
problem would be to solve the stochastic master equa¬ 
tion (SME) numerically for a large ensemble of repe¬ 
titions, and then to perform statistical analysis on the 
sub-ensemble of trajectories which end in li/zf) at time 
T. An alternative approach based on an action principle 
for continuous quantum measurement was developed in 
reference IBl. The action principle naturally incorpo¬ 
rates post-selection and yields a set of ordinary (non¬ 
stochastic) differential equations for the most probable 
path which are simper to solve numerically than the full 
SME. In the limit of no post-selection, this approach is 
consistent with SME formulation. The results of refer¬ 
ence 1^ experimentally verify the predictions of this 
action principle for the case of a single qubit under si¬ 
multaneous measurement and Rabi drive. However, the 
action principle is a general theory which can be applied 
to a wide variety of quantum systems and may prove 
useful in designing optimal quantum control protocols. 


6. Time-symmetric state estimation 

We have so far focused on the use of the quantum 
state as a predictive tool; the quantum trajectories pre¬ 
sented in the previous sections describe the evolution 
of expected average outcome of observables (t_j, ct,,, cTj., 
which relate in a straightforward way to the probability 
of obtaining a certain outcome in subsequent projective 
tomography measurements. However, it is also possible 
to follow the qubit state evolution backward in time to 
predict an unknown measurement result from the past. 



Figure 5: Greyscale histograms of quantum trajectories based on 
5 X 10^ simulated trajectories. Here, t = 1.28jus, y = 2.7 X 10“’ 
and n/2n' = 0.4 MHz. Histograms are normalized such that the most 
frequent value at each time point is 1. Panels A and B depict the full 
distribution of x (A) and y (B) trajectories. Panels C and D display the 
sub-ensemble of trajectories which end in the final state xp = 0.1±.08, 
Zf — .55 ± .08. 


Consider the following guessing game: Two experi¬ 
menters can perform measurements on the same quan¬ 
tum system. At a time t the first experimenter makes 
a measurement of some observable Q.,„ and hides the 
result. The second experimenter then must guess the 
outcome m that the first experimenter received. If the 
second experimenter only has access to the quantum 
system’s state before the first experimenter’s measure¬ 
ment, then the theory of POVMs provides the second 
experimenter with probability for each outcome m and 
the ability to make the best possible guess. However, if 
the second experimenter is allowed to probe the quan¬ 
tum system after the first experimenter has conducted 
her measurements, can he make a better prediction for 
the hidden result? Indeed, since more information about 
the system is available at a later time the second exper¬ 
imenter can make more confident predictions. It can be 
shown MM that the probability of the outcome m is 
given by. 


Pp(m) = 


Tr(Q,nPfQ,'„£'f) 


( 11 ) 


where p, is the system density matrix at time f, condi¬ 
tioned on previous measurement outcomes and propa¬ 
gated forward in time until time f, while E, is a ma- 
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trix which is propagated backwards in time in a similar 
manner and accounts for the time evolution and mea¬ 
surements obtained after time t. The subscript p de¬ 
notes “past”, and in 1461 it was proposed that, if t is in 
the past, the pair of matrices (p,, E,), rather than only 
Pt, is the appropriate object to associate with the state 
of a quantum system at time t. It is worth noting that 
information from before the first experiment’s measure¬ 
ment, which is encoded in p,, and information after this 
measurement (encoded in E,) play a formally equivalent 
role in the prediction for the experimenter’s result. It 
is natural that full measurement records would contain 
more information about the system and several preci¬ 
sion probing theories l47l l48l l49l l50l BTl have incorpo¬ 
rated full measurement records. 

Recent experiments have applied Eq. to systems 
with Rydberg atoms l5^ and superconducting qubits 
Il28ll . confirming how full measurement records al¬ 
low more confident predictions for measurements per¬ 
formed in the past. This applies to both projective and 
weak (weak value) measurements and in the case of 
weak measurements the experiments reveal how the or¬ 
thogonality of initial (p,) and final (£,) states leads to 
the occurrence of anomalous weak values 


7. Two-qubit trajectories 

To this point, we have focused on the trajectories of 
a single-qubit system. However, the trajectory formal¬ 
ism is readily extensible for studying the dynamics of 
a multiple-qubit system. Such systems are of interest 
because they allow us to directly observe and study the 
generation of entanglement on a single-shot basis ll5^ . 
In particular, we study a cascaded CQED system com¬ 
prised of two superconducting qubits housed in two sep¬ 
arate cavities that sequentially probed in reflection by a 
single coherent state (Eigurej^. The quantum trajectory 
formalism for such a system was first developed in 1993 
m, and was observed experimentally in 2014 Il27l . 

The Hamiltonian for this cascaded quantum system is 
given by 


H - Ho+xici^acrl + xib''bcrl - i ^ {a^b-b^ a), 

( 12 ) 

where Hq comprises the uncoupled Hamiltonians for the 
two cavities and qubits; are the dispersive shifts; a (b) 
and ( 1 ?* ) are the creation and annihilation operators 
for the first (second) cavity; /c, are the decay rates of 
the cavities; and rjioss ~ 0.8 represents the transmission 
efficiency between the two cavities. 


In the multiple-qubit regime, it is more convenient 
to work in the measurement basis (|00), |01), |10) and 
111)) rather than the Pauli basis set (cr,- ® crj), since the 
measurement in the multi-qubit case does not project 
along a single-qubit Pauli operator. In the experiments 
described in reference IIZTl . the sequential measurement 
realizes a half-parity operation (Pig 5. A,B): the mea¬ 
surement tone acquires distinct phase shifts of +A0 for 
the even-parity states |00) and |11), and an identical 
(null) phase shift in the odd parity subspace (|01) and 
| 10 ». 

In cascaded quantum systems, the effect of the losses 
between the systems is of primary importance and need 
to be fully taken into account when a quantitative de¬ 
scription is needed Il27l . However in this review, we 
make the choice to set rjioss - 1 (zero losses) for the 
sake of simplicity. In addition, we make the assumption 
that the dispersive shifts are equals for the two cavities 
(Xi - X 2 - z) well as the decay rates (a:i - K 2 - k). 

Similarly to the single-qubit case, we can define the 
measurement outcome 14 as the time-average of the X 2 
quadrature voltage; 14 = IKkAt) Vx 2 (t)dt. The 
dimensionless measurement outcome is thus given by 
r*. = 214 /Ay where AV oc A6 is defined as the distance 
between the measured Gaussian histogram centers for 
|00) and |01) / |10). The measurement realises a projec¬ 
tion on a timescale r s AAtjS with the dimensionless 
measurement strength S - 64x^nri„,At/K. 

The formalism for generating a joint qubit trajectory 
is quite similar to that of the single qubit case. We col¬ 
lect a series of measurements {r<.) at times { 4 ), and use 
these measurements to calculate the conditional density 
matrices {p‘j^'‘"‘], where [ijjm] index the computational 
states. The diagonal density matrix elements can be cal¬ 
culated using a Bayes’ rule, for example: 


^00,00 p00,00|,Q^^[-(rj+2)42<r2] 


(13) 


Here, cr is the width of the Gaussian histograms, 
which decreases as 1/ V^. The off-diagonal density 
matrix elements can also be calculated within the same 
Bayesian formalism. Neglecting internal losses in the 
cavity and Ti relaxation, the off-diagonal density ma¬ 
trix terms p°*'*° are given by: 


Ip 


ijJm. 

k I 


= Ip 


ijJm, 

0 I 


V ijjj lm,lm 

Pk Pk 

V ijjj lm,lm 

Po Po 




X exp 


kAt 

T 




7m 1 2 


(14) 
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Figure 6: Two-qubit trajectories (A) Schematic of the measurement setup. (B) The cascaded cavities are probed with a coherent microwave 
tone, initially aligned along the X[ quadrature. After probing both cavities, the tone can acquire three different phase shifts: +A0, -A6 and 0 
corresponding to the states |11), |00) and |01) / |10) respectively. (C) Conditional quantum state after a measurement result r of length tk = kAt = 
0.8/i5, for both quits initially prepared in the maximally superposed state ((|0> + |1>)/ V(2) <S> (|0> + |1>)/ V(2) and t = 0.15^s. (D) Single quantum 
trajectory of the cascaded two-qubit system. The color code is similar to the one used in panel C. The inset shows the corresponding concurrence. 


where is the average voltage corresponding to the 
qubits prepared in the states i,]. Here, the first term 
represents the Bayesian update and includes intrinsic 
dephasing of the matrix element jijjm', the second 
term accounts for partial dephasing due to uncollected 
measurement photons. Notice that in the case where 
there is no dephasing of the pO'-'o 
diagonal term due to nonunity 77 ,the odd-parity sub¬ 
space becomes protected with respect to added noise, 
and we expect the measurement to probabilistically gen¬ 
erate entanglement in the odd-parity subspace. 

Equations and provide a mapping {r^) 

[pk] at each measurement time f*. = kh.t. We can 
also reconstruct the trajectories experimentally using 
conditional tomography, as described in Section]^ and 
elaborated in reference E?]. Figure 4C shows such a 
Baeysian mapping. However, as mentioned earlier, the 
losses between the cavities need to be accounted for. 
Thus we used a more refined formalism as explained in 
reference IIZTlI . From this mapping, one can reconstruct 
the trajectory of a single iteration of the experiment. 

As mentioned before, the main advantage of this cas¬ 
caded system is its ability to generate entanglement be¬ 
tween remote qubits. To quantify this entanglement, we 


can reconstruct the quantum trajectory of the concur¬ 
rence, which is a monotone of entanglement. A simpli¬ 
fied definition is given by: 

_ o »'v^ov/'r\ / 00,00 11,11 \ /■ic\ 

Ck ~ 2ma.x(0, \p^ \ — yjPj^ Pj^ ) (15) 

Concurrence reaches a maximum value of 1 for a 
maximally-entangled Bell state, and is zero for joint 
qubit states that cannot be distinguished from separa¬ 
ble or from classically mixed states. An exemplar quan¬ 
tum trajectory of the joint qubit state, and of the con¬ 
currence, are shown in Figure 4D. We see an initial 
transient during which C is zero, followed by a non¬ 
monotonic increase in the C as the joint qubit state 
stochastically projects towards the entangled manifold, 
reaching an eventual concurrence of C ~ 0.55, indicat¬ 
ing a highly nonclassical state. 

8. Outlook 

The experiments presented in this review demonstrate 
precise control and a detailed understanding of the pro¬ 
cess of continuous quantum measurement of a super¬ 
conducting qubit. This knowledge may benefit a wide 
range of future research directions in quantum control 
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and multi-qubit state estimation ll54l l55l . In this fi¬ 
nal section we highlight one such research direction: 
measurement-based quantum feedback. 

A continuous measurement record, which contains 
information about how the qubit state evolves in real 
time, can be incorporated into a feedback loop for a 
number of applications including state preparation, state 
stabilization, and continuous quantum error correction. 
Without feedback, a combination of projective measure¬ 
ment and unitary rotation can be used to probabilisti¬ 
cally prepare an arbitrary qubit state ll5^ . Using feed¬ 
back, the measurement result can be used to control 
a subsequent qubit rotation, allowing for deterministic 
state-reset protocols ll57l |58l which can be repeated on 
a timescale much faster than T i. 

In the case of weak measurement, it is possible pre¬ 
pare arbitrary states by measurement alone, without ap¬ 
plying any subsequent qubit rotations. Without feed¬ 
back, state preparation is probabilistic: one simply post¬ 
selects an ensemble of trajectories which end up in the 
desired state. With feedback, it is possible to prepare an 
arbitrary qubit state deterministically through adaptive 
measurement, as recently demonstrated using nitrogen 
vacancy centers |[33l. 

In addition to state preparation, quantum feedback 
can also be used to stabilize a qubit state or trajectory. 
In reference 1301 . we used weak measurements and con¬ 
tinuous quantum feedback to stabilize Rabi oscillations 
of a superconducting qubit. Reference ll58l demon¬ 
strated that stroboscopic projective measurements and 
feedback can be used to stabilize an arbitrary trajec¬ 
tory, such as Rabi or Ramsey oscillations. Furthermore, 
with the with the two-qubit setup from reference Ezl 
it should be possible to use feedback to to stabilize an 
entangled state. 

Looking forward, proposals for fault tolerant quan¬ 
tum rely on quantum error correction (QEC) protocols 
in which a single logical qubit is composed of many 
physical qubits. While many QEC schemes such as sur¬ 
face codes rely on discrete projective measurements of 
syndrome qubits iB^lbOllbTl . a wide body of QEC pro¬ 
posals are based instead on continuous measurement- 
based quantum feedback In these techniques, a sin¬ 
gle logical qubit is encoded in several physical qubits, 
and an error syndrome is detected by processing (or ‘fil¬ 
tering’) a continuous measurement signal. The error 
signal is used to generate a suitable feedback Hamil¬ 
tonian which corrects for errors in real time. 

By tomographically validating individual quantum 
trajectories, the experiments presented in this review 
have demonstrated the ability to correctly ‘filter’ a mea¬ 
surement signal for one and two qubit systems. A sen¬ 


sible next step is to build a system of several qubits and 
attempt to correctly filter an error syndrome. Then, the 
following step would be to feed-back on this error syn¬ 
drome to realize a single logical qubit whose lifetime 
exceeds that of its constituent physical qubits. To re¬ 
alize this goal will require a robust multi-qubit archi¬ 
tecture and improvements in the measurement quantum 
efficiency. Recent experiments have demon¬ 

strated that it is possible to individually measure and 
control 5-9 qubits in a planar cQED architecture, and 
efforts to improve the measurement quantum efficiency 
are currently underway in a number of different research 
groups. Although there are still formidable challenges 
to overcome, and while the ultimate utility of measure¬ 
ment based QEC in comparison to other methods re¬ 
mains an open question, it seems that an initial demon¬ 
stration of measurement based QEC may lie on the hori¬ 
zon. 
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